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ABSTRACT 

A  pertubation  expansion  is  formulated  for  the  three- 
dimensional,  nonlinear,  acoustic-wave  equation  with  dissipa- 
tive  term  describing  the  viscous  and  thermal  energy  losses 
encountered  in  a  cylindrical  cavity.   The  theoretical  results 
show  that  nonlinear  effects  in  sloshing  modes  are  strongly 
suppressed . 
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LIST    OF    SYMBOLS 


c 

Co 

Ci 


Po 


u 

If 
id 


phase  velocity  in  the  cavity 

=  cV  -  a'/at* 

effective  phase  speed  associated  with  a  standing  wave 

at  resonance 

specific  heat  at  constant  pressure 

Wy£,  propagation  constant  associated  with  a  standing 

wave 

peak  Mach  number  of  the  driven  standing  wave 


P 


-  \J    -    -  \o   ^-/^f  =  acoustic  ores  sure 


JL 


instantaneous  and  equilibrium  total  pressure  in  the 

field 

absolute  temperature  in  kelvin 

particle  velocity 

particle  speed  of  the  Tl^h-order  perturbation  solution 

infinitesimal-amplitude  attenuation  constant 

temporal  decay  constant  of  a  resonance 

=(T  +  1)  /2  for  a  gas 

=  ^-P/Cv5  ratio  of  specific  heats  for  a  gas 

instantaneous  and  equilibrium  densities 

velocity  potential 

=aXt,  (angular)  grequency  at  which  the  cavity  is  driven 

(angular)  frequency  of  a  resonance 


I.   INTRODUCTION 

The  topic  of  finite  amplitude  acoustic  standing  waves  in 
sloshing  modes  of  a  cylindrical  cavity  is  interesting 
theoretically,  but  development  of  the  subject  has  not  been 
extensive.   The  purpose  of  this  research  is  to  present  the 
results  of  a  power  series  perturbation  approach  to  the 
problem. 


II.   THE  NON- LINEAR  WAVE  EQUATION 

A.   GENERAL 

It  is  well  known  [l]  that  for  °^/\ and  M«l,  where 
measures   the  fractional  loss  per  wavelength  and  M  is  the  peak 
Mach  number  of  the  source,  loss  terms  and  nonlinear  terms  in 
the  constitutive  equations  can  be  separately  approximated  with 
the  help  of  linear,  lossless  acoustic  relations.   The  force 
equation  appropriate  for  acoustical  processes  in  systems  for 
which  gravitational  effects  are  unimportant  is 


at 


c  u-v)  u+  e"Vp  =^*cu  (2-D 


where  (_J  =  Li  [^\   is  the  particle  velocity,  pis  the  instantaneous 
density  of  the  fluid,  P  is  the  instantaneous  total  pressure  in 
the  field,  and  the  operator  ^symbolically  describes  those 
physical  processes  leading  to  absorption  and  dispersion.   We 
used  two  additional  equations. ■  The  first  is  the  equation  of 
state  for  a  perfect  gas 

P   ~   (T~k  (2-l-a) 

where  r  is  a  constant  whose  value  depends  on  the  particular 
gas  involved  and   |K  is  the  absolute  temperature  in  Kelvin. 
The  second  is  the  continuity  equation 


4-  +  v-[<  i+s)u] 


-   O 


(2-1-b) 


where  s  =  (  P~P«>)  /  ?°     is  the  condensation  at  any  point  and  P 
is  the  equilibrium  density  of  the  fluid.   If  we  ignore 


rotational  effects,  then 


u    =71 


(2-2) 


where  O  is  the  velocity  potential.   Combination  of  Eqs.(2-D- 
(2-2)  and  the  neglect  of  terms  of  orders  higher  then  M  ,  M(°V-j^)  , 
and  (   A&)  ,  yields  a  quadratically  nonlinear  wave  equation, 

P   * 
3t  ^'^G* 


(2-3) 


-^ 


D. 


where  Co^  =  (  i/^p  )  (adiabatic),  p  =  acoustic  pressure ,M  =(J  .  M 
^  is  the  ratio  of  heat  capacities,  and 


c;n*  =  c:v*- 


^ 


(2-4) 


The  left-hand  side  of  Eq.(2-3)  is  the  classical,  linear  wave 
equation  with  losses  pertinent  to  the  system  under  study. 


The  right-hand  side  can  be  interpreted  as  a  forcing  function 
consisting  of  a  three-dimensional  spatial  distribution  of  phase- 
coherent  sources. 

In  a  second-order  perturbation  theory,  this  volume  forcing 
function  is   obtained  from  the  classical ( first-order )  solution 
P]_  of  the  acoustical  problem.   The  second-order  perturbation 
solution  P]_  +  P2  describes  the  non-linearities  P2  resulting  from 
the  self  interaction  of  the  classical  solution  P-j_. 
Higher-order  perturbation  solutions  consider  the  interaction 
of  the  nonlinear  solution  with  itself,  and  the  forcing  function 
is  composed  of  products  of  both  classical  and  non-linearly 
generated  terms.   Thus,  if  a  system  is  driven  at  frequency  w, 
the  non-linear  term  in  Eq.(2-3)  will  force  the  existence  of 
all  integer  multiples  T)U/  of  the  driving  frequeney  and  the  full 
solution  must  contain  all   harmonics  of  the  input  frequency. 
In  a  closed  cavity,  each  of  those  nonlinearly   generated  waves 
whose  frequency  lies  near  the  resonance  frequency  of  a  standing 
wave  of  the  cavity  and  whose  assciated  spatial  function  matches 
that  of  the  standing  wave  can  be  strongly  excited  [2]  . 

As  far  as  the  author  has  been  able  to  determine,  there  has 
been  only  one  previous  study  of  this  system  published  in  the 
open  literature.   This  was  by  Maslen  and  Moore  in  19  5  6  [3]  . 
Their  approach  resulted  in  a  series  expansion  which  did  not 
converge  if  the  relevant  nomal  mode  frequencies  were  integerally 
related.   Their  interpretation  of  the  quenching  of  the  nonlinear 
effect  was  based  on  the  "scattering  effect  of  the  wall". 


We  feel  our  interpretation  based  upon  nonlinearally  generated 
volume  sources  stimulating  the  allowed  standing  waves  of  the 
cavity  is  more  accessible  and  informative.   Further,  the 
mathematical  approach  developed  herein  appear  to  avoid  any 
difficulties  in  convergence.   Their  conclusion   that  high 
amplitude  monofrequency  transverse  oscillations  can  exist  is 
consistent  with  our  finding. 


B.   APPLICATION  TO  THE  CYLINDRICAL  CAVITY 

The  circular  cylinder  is  provided  with  a  "point"  source  of 
sound.   By  properly  positioning  the  cylinder  with  respect  to 
the  sound  source  it  is  possible  to  effectively  drive  the 
enclosed  air  into  various  modes  of  vibration.   The  rigid- 
walled  cylinder  forces  the  component  of  particle  velocity 
perpendicular  to  each  cavity  surface  to  vanish  at  the  surface. 
The  resulting  steady-state  solution  to  the  linear  wave  equation 
in  cylindrical  coordinates  is 


D 


(  T,  $ .  J )  =  Amni0*  ^  *  ^ m  ^V006  OW^ InC*nn*r) 


(2-5) 


where  Jn  are  the  cylindrical  Bessel  functions  and  application 
of  the  boundary  condition  to  the  sides  yields 


4 


_    J 


WW) 


7i -m 
CL 


(2-5-a) 


10 


where  a  is  the  radius  of  cylinder  and  J^  are  the  arguments 

ik 

of  the  exterema  of  the  7]   Bessel  function. 

The  nomal  mode  frequency  is  dependent  on  -Kay  and  "KyrVW  -> 


^ 


X. 


o     _«_  vvjTvnyL 


\A, 


*««,    +  %      -       G*  (2"5-b) 

The  standing  wave  well  be  identified  by  the  ordered  integers 
(n,  m,  X  )  describing  its  spatial  dependence. 
1 .   Symmetric  Modes 

The  simplest  waves  in  cylindrical  coordinates  are  those 
that  depend  only  on  the  distance  r  from  Z-axis,  the  gradient 
takes  the  form 

y=$~k±z  (2_6) 


A 


where 


^  =  -kr 


(2-6-a) 


and  k=w/c  is  the  wave  number  or  propagation  constant 
The  Laplacian  becomes 


a?        i   3 


XJ     -  ~k    (  Mj*  +  *j  ?|  )  (2-6-b) 

The  suitable  trial  solution  appropriate  for  symmetric  modes  is 
the  linear  solutionCO,  m,  0)  of  frequency  w 


11 


ll  =  m  X  ^)  /m-  ^ 


(2-7) 


w 


The  velocity  potential  for  this  limit  can  be  defined  as 


Co' 


1  =  C  ^  )  ^  Jo  C^)  cos  ukt 


(2-8) 


,11,, 

and         I r      is 

u  _  " 


G 


i  M  Jo^)  cos  tot  =- |M  J,cc6  i*rt 


(2-9) 


Because  we  assume  that  the  surfaces  of  the  cavity  are  rigid, 
then  at  r=a  the  appropriate  boundary  condition  is  Jq'=0. 
Thus,  we  have  ~%(k.   =  J  q^  • 
To  generate  the  second  order  solution,  we  first  note  that 


cM-±tfj.zo -«*»*) 


(2-10) 


and 


Ll 


Z-r-Z 


r  i^_  y=  ±  p]  J(r  ( l  +  cc6  zurt) 


(2-11) 


Next,  the  second  derivatives  of  Eqs.(2-10)  and  (2-11)  with 
respect  to  time  are 


12 


C—  Mlr«s«rt 


(2-12) 


and 


2_ 


U,    x-  (*W    N/i^-r'1 


(  ~)   = 


G 


/\^ 


|vj  j  (  cos  zjtft 


(2-13) 


Substituting  Eqs.(2-2)  and  (2-13)  into  the  first  term  of  the 
right-hand  side  of  Eq.(2-3)  yields 


ssf  m-cos  ™t\  r-  rj:} 


(2-14) 


and  the  second  term  of  right-hand  side  of  £q.(2-3)  becomes 


(2-15) 


since 


Z         -t-Z 


c:/r-^(3r-J. ) 


(2-16) 


or 


Co  V  J,  =*w  CtX  +  *  J.  -J,  ) 


(2-17) 


With  the  use  of  Eqs.(2-14)  and  (2-15)  we  can  write  the  right- 


13 


hand    side    of    Eq.(2-3)    as 


4 


>        f 


cos  Mot  y 


t 


■*  .  ~rz     rv  -L 


1I  i-^  r-r-l-rj  +  -ir+r-a-t)j. 


• 


(2-18) 


and   the  appropriate  inhomogeneous  wave  equation  for  the  second 
order  perturbation  solution  is 


Re  (2W)J 


i'      I  —  z 


ccdHeJC^  -  ^■ri'[[-iJ."+J,-tr]+(-iI 


H- 


J-Cr-^j/jcesxwil 


(2-19) 


2 .   Non- Symmetric  Modes 

The  non  planar  waves  in  cylindrical  coordinates  are 
those  that  depend  on  the  distance  r  from  Z  axis  and  the  angle 
Q   from  X-axis.   The  gradient  takes  the  form 


v-=^-ci!^i# 


(2-20) 


where 


1=    ^r 
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and   the    Laplacian   becomes 


i 


V*      =   "fe     ^If    +y  9$    +     f  Py  ^  (2-20-a) 

The  suitable  trial  solution  appropriate  to  the  forcing  function 
of  frequency  w  for  exciting  the  (n,  m,  0)  standing  wave  is 

-H_  =  m  j^tf^r^namrtt  r  -few- j/iw/a.  (2-21) 

where  Jn(xD  is  the  Bessel  function  of  the  first  kind  and 
order  n.   The  velocity  potential  is  approximated  by 

J  =  (l£-)  [""J  X,  C-fer)  0*7)$  06  lift  (2-22) 


ind  U«  /  /*   is 

U*_  -  [V](  ^  J^CoS^COSWft-^  ^X^^CCStftj  (2-23) 
Co 


Thus, 


Hjg)  =  M^jJ^5^  ^^  (2_24) 


or 


f  II  f  =-^MxC  \-<&™tXX+T*Xcc*  Z^V        (2"25) 
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and 


rW-i  M2d  +eosait)(i(^+ j;  l-VL^  x^j 


G^    4 


(2-26) 


and  the  second  derivative  of  Eqs.(2-25)  and  (2-26)  with  respect 
to  time  is 


■f   /  _P^  <*_  &M  M  *  J*  (  |  +60S  iDfi  )  CCS  XlM-t 


tf^'      4 


and 


1 ,  u, 


3t 


Uo 


=.^(^+(1)3 


7)  nz  -r*l 

71 


+ 


r  -w*  ('7i\aTaN 


T'  _  j-lL)  I    COS  */)!$ 


(2-27) 


(2-28) 


Substituting  Eqs.(2-26)  and  (2-27)  into  the  first  term  of 
Eq. ( 2-3 )  yields 


(2-29) 
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or 


T*l 


(2-30) 


since 


j^  =  ^Cj^  +  J^,^T^h) 


J 


(2-31) 


and 


on   T* 


L/  lC-^3^  +  xJwi^w 


+ 


(2-32) 


In  order  to  solve  the  second  term  of  Eq.(2-3),  we  can  use 

z^z     -,  *.*/     '    T^  IT*  -T  )  (2-33) 


(2-34) 


(2-35) 
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since 


J 


x. 


(  T  )  J-^  =  T  (  3tw  '     Jm  ) 


3 


(2-36) 


and 


J*"  (-S-)     J^     -         JjhH  J^-i 


(2-37) 


Thus,  Substituting  Eqs.(2-33)  through ( 2-35 )  into  the  second 
term  of   Eq.(2-3)  yields 


f 


c*V'(-^)*=+M,lvtfiCl-«6xWrt)[c43 


# 


3  >. 

JiH  ~X  )  +  (-  3™  J*h  ~  X  )  COS  am^ 


(2-38) 


and 


J--T 


-  3^)  ~  T41+1  *+"  z,  Jin*     J^i  /  ^v  z,  J 


+■  £  1***3**  +  J™  Xi  )c^  ^ 


(2-39) 


With  Eqs.C2-30),  (2-38)  and  (2-39)  we  can  get  the  inhomogeneous 
wave  equation  for  the  second  order  perturbation  as 


c^n^ftx  )!==  fV[fu:+x:-ir 


$G"  4 


4  Jin*        4  3*-,*,  )  T  4  (  ZJttiJt- 


N 


.1  / 


z 


(M) 


4_  ^  ri  ces  ^i*t    [-jj3i«  +  ^Jm 


4- 


-(X.^+J-^ 


4.  ^H         |6 


■J*, 


X 


-  g-)j;Jcos^ 


(2-40) 
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If  we  let  n=0  in  Eq(2-40),  it  reduces  to  Eq.(2-19),  as  it  must 
Substituting  n  =  l  into  Eq.(2-Li-0)  yields  the  equation  with 
forcing  term  resulting  from  the  (1,  m,  0)  sloshing  mode, 


p.         <>) 


C<^X+-£jC)#-  ^  M ^T(-?Tf+J. 


^-r;+fc-lxv^j3j 


> 


*J*Jo  )<»6  ^J  +^  M^^wH:. 


J 


r 


6  J3 


4- 


X       (  x  i/  )  J;    +4.3o 


*>  |6 


f 


.> 


t^^Ti   "  TJ-I~(~J""^)3f   GCfc^] 


> 


(2-41) 
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III.   METHOD  OF  SOLUTION 


Recall  that  the  equation  with  forcing  term  resulting  from 
the(l,  m,  0)  sloshing  mode  can  be  written  as  a  function  of 
frequency , 


4 


\L 


Ta4-J-T*-P^-J-n+-L-r 


-t 


V. 


TZ 


{lu-jjjc-^-T'Jrj^^j 


(3-1) 


and  the  left-hand  side  of  Ea.(2-3)  is 


\   r  z  jZ  /  J> ,  J — L 


S 


P. 


?i   '  ?t 


'J  Po^o 


(3-1-a) 


If  the  harmonics  of  the  frequency  at  which  the  cavity  is  driven 
are  not  close  to  any  resonant  frequencies,  we  can  ignore  the 
lossy  term.   So,  Eq. (3-1-a)  can  be  written  as 


*.*/  31 


I   9 


&2"fc  C*-n*+  ^  H  ) 


n  t  nJ  *us 


i 


(3-2) 


i  Co 


Let  us  assume  that  the  solution  of  Eq.(3-1)  can  be  expressed  as 
P 


^oCo 


s  =    rfccS  Xtftl    UsntyHM^V  "*■  X75 


(3-3) 
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and   define 


P 


e 


311  =  f|   V^  Ofay)  cos  -zrisQ  ccGziut 


iUJ 


(3-4) 


and 


P,Cc" 


6 


(3-5) 


Combination  of  Eqs.C3-2)  and  (3-3)  yields 


\ 


^ 


LUc^)= 


/ 


4 


-  — T    4— Hi 


jr 


s 


o      J«  +"     4-  Jtjh       \l  J- 


f  Jth     /fcJoi-^j 


(3-6) 


and 


\    XW  J 


t\ 


^+^-(T)J+,r^-4ttw 


^T    T    _GL__Ut 


J. 


(3-7) 
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A.   POWER  SERIES  METHOD 

In  order  to  solve  the  right-hand  side  of  Eq.(3-7),  we  can 
use  the  definition   [  4 


1)  ^~ 


Jv^)XtCf  =  l^ 


J,=o 


*M 


(3-8) 


where    [""*  (  J-j-  I  )  =  ^  J 

1.   Coefficient 


a.   Substituting  V  =  ^  =  n+2  into  Eq.(3-8)  yields 


J«^)=(^) 


W4.^« 


■fcO 


Hf  (^*fc^)/("zi) 


4? 


>  *■ 


(3-9) 


or 


A3 

T^(1)=X    Brvfeb 


i(fH--feivt) 


■4=0 


(3-10) 


where 


B**=(-< 


4g   (^^z-fet  4)  / 


(iI-Mm-*)  /       C*rt+-fet4.)/-fe/ 


(3-11) 


b.       Substituting    V=    n+2,^  =  n-2    into    Eq.(3-8)    yields 


vWC^)J^,^)  =  C— •  i, 


ypTft 


fin 


2l 


-H 


** 


HfC^^)/(^1) 


2* 


(1I-+4+*)/  (^t4!-x)/(^^)/^/ 


(3-12) 


or 


I^(l)T^) 


Co  i^) 


4=0 


(3-13) 


where 


(3-11) 


c.       Substituting  V)  =  n+1, .{(  =n-l    into    Eq.(3-8)    yields 


Jw(1)Ih(1)=(*1 


iH 


CX> 


-frO 


(3-15) 


or 


■fro 


X^f-fe) 


(3-16) 
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where 


where 


(3-17) 


d.   Substituting  V=  {! =n  into  Eq.(3-8)  yields 


(3-18) 


or 


Co 


■fro 


i(OH-fe) 


(3-19) 


t/ft£ 


(3-20) 


e.   Substituting  V=^( =n+l  into  Eq.(3-8)  yields 


x+l(^=Oi 


i^l+Z.  <>° 


■* 


M-ft 


:if(zwz|tx)l  (-^J 


>* 


>x 


(nH-Hl)!  (x-n-r^O  !  * 


(3-21) 


25 


or 


c*=> 


.c-nt^ti; 


(3-22) 


where 


^Tt*   ^ 


[(TVt-fc+l)/   C^l+M)!*: 


(3-23) 


f.   Substituting  V  =^(.  =  n-l  into  Eq.(3-8)  yields 


or 


XVl)  =  (iD" 


Zf-Z 


<>=> 


*o 


p  lie 


Co 


4=0 


iC'TH-'fe-f ) 


(3-24) 


(3-25) 


where 


H^  =e^ 


(  a^]^a^-^) ! 


[<M^)\]\*^olM 


(3-26) 


26 


g.   Substituting  ])  =/ii  =n-2    into  Eq.(3-8)  yields 


or 


#0 


Oo 


-ho 


Ual£A*1J 


2. 


(th-fc-zO 


(3-27) 


(3-28) 


where 


.i 


out-t-'fr.  ,fjr,V 


(111-fc-zO  !   [  (*m>fe-4)/-fe/ 


(3-29) 


2.       Vn(y) 


With  the  help  of  Eqs.(3-9)  through ( 3-2  9)  Eq.(3-7)  can 
be  written  as 


XA 


iL.+  _L_L 


£2L) 

1/ 


z^ 


+ 


V»c^) 


oo 


(3-30) 


27 


where 


x 

Z 


(3-30-a) 


and 


a  = 


(3-30-b) 


Now,  let  us  assume  that 


-fe) 


(3-31) 


The  first  derivative  with  respect  to  y=-fer  can  be  expressed  as 


'1 


1- 


y~  i^-k)  cucii 


:i(Hl~fe)-' 


(3-32) 


Hh 


and  the  second  derivative  with  respect  to  y=^r  can  be  expressed 


as 


d*V-  _ 


Ji 


i 


,i(^)-X 


(3-33) 


Substituting  Eqs.(3-31)  through( 3-33 )  into  the  left-hand  side 


of    Eq. (3-30)    yields 


co 


^ 


( »)*  iC^^K^-fc-i)  CUc-^5 


X<7)+4H ) 


?0 


28 


4- 


£e>i++)  GUC^l) 


2(<rH-£)-i 


X^t-JO' 


-muM-Hf  j+w*ij 


jjit-k) 


(3-34) 


Eq.(3-34)  must  be  equal  the  right-hand  side  of  Eq.(3-30).  Thus, 


aa 


)Qmk+ 1  X-fet-'Xl  hi  +GU=  ~  ("8^* + i  D*"1"  A  E 


(3-35) 


where    k    =    0,    1,    2,    3       (3-36) 

If  Eq.(3-35)  is  substituted  into  Eq.(3-36),  then  Eq.(3-35)  can 
be  expressed  as 


a«  +  (irj C ^+3)(  * )  CU=-(iCtt+iD»-f  AE«t) 


(kt  C^ 


(±T  c^-fe+i  )(^+| )  ^>=  ~(K*+  ^  ^+  ^  £•*< 


(3-37) 
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The  normal  component  of  particle  velocity  is  equal  to  zero  on 
the  boundaries.   From  the  application  of  the  boundary  condition 
at  the  rigid  boundary  at  r=a, 


<JV« 


d  ^ 


s-o 


i 


(i»+*)(iU(»j' 


aClH--fe.>l 


=  0 


IT 


UrO 


where 


JC-K)  =  <ll/<A(-K> 

Substituting  Eq.(3-36)  into  Eq.(3-38)  yields 


(3-38) 


-K*n*)OU  (ijO^  +  +(^)(WV3 


,-^6w4h) 


=  O 


(3-39) 


3.   Un(y) 


With  Eqs.(3-9)  through( 3-29 )  ,  Eq.(3-6)  can  be  expressed  as 


-fcO 


AtfH-fe) 


; 


X-r+fe) 


■+"  4    Hh,*h  (*|)       ~li©^(^ 


-2(<7H4-Z)N 


(3-40) 


30 


where,  as.  before,  A  =  0/2-1/8  and  a  =  T"j£ 
Now  let  us  assume  that 


4 


(3-41) 


^=0 


The  first  derivative  with  respect  to  y=-fcr  can  be  expressed  as 

(3-42) 


and  the  second  derivative  with  respect  to  y=-j^frcan  be  expressed 


as 


lit  _  J.  V  ciH*4Kwt^)UkC^- 


(3-43) 


Substituting  Eqs.(3-41)  through( 3-4 3 )  into  the  left-hand  side  of 
Eq.(3-40)  yields 


f 


-fee 


(il+40  ~4<^-^H ) 


s 


^-fc  \  ^  c 


z1 


"J-  L-&.  (  *0) 


.CTH-^-Z)  / 


(3-44) 


Eq.(3-44)  must  equal  the  right-hand  side  of  Eq.(3-40). 
Thus, 


31 


{lk1 'OiH-k-t-S  ?"■■&>  +  •*!*  =  (—}£&».*+  +f^.*-> 


fU.-i&O-AE 


4.  1  iTJ.-fc-l       16 


?.•&-> 


-/v 


(3-45) 


where  ^  =    1,    2,    3         

If  Eq.(3-45)  is  substituted  into  Eq.(3-36),  then  Eq.(3-45)  can 
be  expressed  as 

w^rd      — La 


k,+(i;)^-^U.    =-■£& 


1™  +(i)1C/"  -  if  Ut  -"  it  &w  +4  H» 


.L/i 


(3-46) 
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Application  of  the  boundary  condition  at  r=a  yields 


where 


Cut*-*)  U  (*J )       =0 


^C-R^J/deH) 

Substituting   Eq.(3-36)    into    Eq.(3-47)    yields 


^JHt) 


p/nc  v.  X 


tC^U^*/)    4  (ohokGj 


+ 


(3-47) 


(3-43) 
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IV.   CONCLUSIONS 
Recall  that  we  can  compute Km^  from  Eqs.(3-37)  through ( 3-39 ) 
and  Eqs.(3-14),  (3-17),  (3-20).   Thus,  with  the  use  of  Eq.(3-31), 
we  can  get  V  .   From  Eqs.(3-46)  through( 3-48 )  and  Eqs.(3-ll), 
(3-20),  (3-2  3),  (3-26),  (3-29),  we  can  compute  [\^   .      Thus,  with 
the  use  of  Eq. (3-14) ,  we  can  get  Un.   Therefore  we  can  compute 
Vn(y)  and  Un(y).   Let  us  calculate  the  finite  amplitude  effects 
resulting  from  the  nonlinear  distortion  of  a  forced  radial  mode . 
If  we  excite  a(0,  m,  0)  mode  and  obtain  the  pressure  at  the 
circumference,  then  p.  ■  =Jfecm&  =  J^  •   Given  this  value  of  M    ,  use 
of  Eqs.(3-14),  (3-17),  (3-20)  give  the  quantities  C,  D,  E. 

A 

From  Eqs.(3-35)  and  (3-36)  we  can  compute  U'flfc  •   Thus,  we  can  get 
Vn(y)  from  Eq.(3-31).   Let  us  calculate  the  finite  amplitude 
effects  resulting  from  the  nonlinear  distortion  of  a  forced 
sloshing  mode.   If  we  excite  a(l,  m,  0)  mode  and  obtain  the  pre- 
ssure at  the  circumference,  then  ^J  =  -y^Q.  -    \y^    .   Given  this  value 
of  4   ,  use  of  Eqs.(3-ll),  (3-20),  (3-23),  (3-26),  (3-29) 


g;ive  the 


quantities  B.  E,  G,  H,  0.   From  Eqs.(3-45)  we  can  compute  j^£ 
Thus,  we  can  get  Un(y)  from  Eq.(3-41).   Substituting  n=0  at 
radial  modes  into  Vn(y)  and  Un(y)  yields  Vg(y)=Ug(y).   Therefore, 
with  the  use  of  Eq.(3-3),  we  can  get  the  solution  of  Eq.(3-2). 
A.   RADIAL  MODES (0,  m,  0) 


34 


then 


=  zhfVoCf^cs  xwt 


so    that 


P>!  l_ 


IP-I 


=  x  ii  Vo  c  yo«,  V  joCj"^) 


r»i 


B.   SLOSHING  MODES (1,  m,  0) 


RCd* 


then 


so  that 


iL=tifV0(jj+ViCJN)H/LJ'cJ- 


p< 


r=& 
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